GW band structure of InAs and GaAs in the wurtzite phase 
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We report the first quasiparticle calculations of the newly observed wurtzite polymorph of InAs 
and GaAs. The calculations are performed in the GW approximation using plane waves and pseu- 
dopotentials. For comparison we also report the study of the zinc-blende phase within the same 
approximations. In the InAs compound the In 4d electrons play a very important role: whether 
they are frozen in the core or not, leads either to a correct or a wrong band ordering (negative 
gap) within the Local Density Appproximation (LDA) . We have calculated the GW band structure 
in both cases. In the first approach, we have estimated the correction to the pd repulsion calcu- 
lated within the LDA and included this effect in the calculation of the GW corrections to the LDA 
spectrum. In the second case, we circumvent the negative gap problem by first using the screened 
exchange approximation and then calculating the GW corrections starting from the so obtained 
eigenvalues and eigenfunctions. This approach leads to a more realistic band-structure and was 
also used for GaAs. For both InAs and GaAs in the wurtzite phase we predict an increase of the 
quasiparticle gap with respect to the zinc-blende polytype. 

PACS numbers: 71.15.Qe, 71.55.Eq, 71.20.-b 
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With the development of the growth techniques such 
as Chemical Beam Epitaxy (CBE) new nanostructured 
materials are being created. Among the technologically 
important III-V semiconductors, one exceptional exam- 
ple is InAs. Indeed InAs nanowires (NWs) grow purely 
in wurtzite structure with [OOOllorientation when InAs 
(lll)B is chosen as substrate [l|. Recalling that the 
zinc-blende structure (zb, 3C, space group F43m (Tj)) is 
the stable phase of bulk InAs, the new hexagonal phase 
(wurtzite wz, 2H, space group P63UIC (Cf^)) clearly rep- 
resents a theoretical challenge. Obviously, the new phase 
calls for a proper theoretical investigation in order to bet- 
ter understand its physical properties and give a correct 
interpretation of recent experiments conducted on InAs- 
based NWs. Some examples are given by the photolumi- 
nescence (PL) measurements performed by one of the au- 
thors on InAs-InP core-shell strained NWs Q , photocur- 
rent measurements on InAs NWs with a centrally placed 
InAsi^Pz segment Q, and electron g-factor measure- 
ments on InAs NW quantum dots [1] . GaAs NWs grown 
on GaAs (lll)B substrate, instead, typically consist of 
alternating wurtzite and zinc-blende segments. 

The electronic structures of the wurtzite polytypes of 
InAs and GaAs are almost unknown. This is certainly 
true for the variation of the fundamental gap of the poly- 
type holds in particular for the varation of the fundamen- 
tal gap with the polytype. Moreover, once the full band 
structure of the new polymorph is obtained, it can be 
used to extract parameters needed for tight binding and 
k-p perturbation theory calculations. Furthermore, the 
so obtained k-p Hamiltonian is a necessary input to elec- 
tron transport studies on the same material, with obvious 
consequences to nano-elcctronics. 



Since there are no measured properties of InAs in 
the wurtzite phase, the choice of an ab initio method 
is highly appropriate because it does not require any 
previous knowledge of the material parameters and al- 
lows one to predict numerous properties ranging from 
the equilibrium lattice constant to the band structure. 
The ground state electronic structure of solids is well de- 
scribed by Density Functional Theory (DFT) Q in the 
Local Density Approximation (LDA) Q. Hence, these 
calculations are specifically suited for structural stud- 
ies. However, the LDA is not appropriate for describ- 
ing electronic excitations. In particular, the calculated 
band gaps are underestimated. The excited-state prop- 
erties of a many-electron system, such as the band struc- 
ture, require a quasiparticle calculation in order to prop- 
erly account for many-body effects. These calculations 
rely on the concept of the single-particle Green func- 
tion G, whose exact determination requires the complete 
knowledge of the electronic self-energy E. The latter is 
a non-Hcrmitian, non-local, and energy-dependent op- 
erator which describes exchange and correlation effects 
beyond the Hartree approximation and which can only 
be calculated approximately. A useful approximation is 
the dynamically screened exchange or GW approxima- 
tion (GWA) 0, @|- Since the electron density is accu- 
rately given by the LDA, the wave functions of this ap- 
proximation are usually rather close to the quasiparticle 
wave functions obtained from the self-energy operator. 
Therefore, the quasiparticle band structure is usually cal- 
culated pcrturbativcly starting from the LDA wave func- 
tions and eigenenergics. 

Our goal is the GW band structure of InAs and GaAs 
in the wurtzite phase. In the case of InAs, the main ob- 
stacle on this path is the fact that the correct procedure 
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to include the In Ad electrons among the valence states 
will erroneously predict InAs to be a metal within the 
LDA. Hence, the LDA cigenfunctions and energies are 
no longer a good starting point for the GW calculation. 
The main reason for the failure of the LDA is an over- 
estimation of the pd coup ling as has been concluded by 
several authors fol. fTol. [Til . Il2f in different materials like, 
for instance, InN [l3| . GaAs [l4| and II- VI compounds 
[l5j . If, instead, the In Ad states are frozen into the core, 
the correct band ordering can be re-established, making 
the LDA results a reasonable starting point for a subse- 
quent GW calculation. For this specific reason, in the 
first part of this study, we have used In pseudopotentials 
for which the d electrons are treated as core states and 
their self-interaction corrections (SICs) are included in 
the underlying free-atom calculation [161 ]. 

Another way of circumventing the negative gap prob- 
lem, is to perform a calculation using an approximate 
self-energy which restores the correct band ordering. The 
screened-exchange approximation is such a self-energy 
The GW calculation is then based on the latter results. 
The prf-repulsion is correctly taken into account. 

Even if there is no problem with the band ordering, 
as is the case in GaAs, it could still be valuable to base 
the GW calculation on an electronic structure obtained 
from a self-energy producing results much closer to the 
GWA. This would lead to final results closer to a self- 
consistent GW calculation. We have found this to be 
important also in GaAs. Previous calculations treating 
the GW self-energy as a perturbation to the LDA poten- 
tial have resulted a too small band gap. This is true also 
in all-electron calculations [13] ■ Thus, it appears that 
in materials with d-electrons like, for instance, InAs and 
GaAs, self-consistency is an important issue. 

Because of the novelty of the wurtzite polymorph, 
there are few theoretical studies at the DFT/LDA level 
fl8L [IH reported in the literature and no GW calculation. 
Besides, it should also be noted that these DFT/LDA cal- 
culations were done assuming the In Ad electrons frozen 
in the core and this, as we will discuss in this article, leads 
to an underestimation of the equilibrium lattice parame- 
ters and to a misleading description of the band-structure 
if the pd repulsion is not properly taken into account. 

On the basis of total-energy calculations within DFT 
it has been concluded [2(| that the wurtzite phase is a 
metastablc high-pressure modification of the InAs and 
GaAs compounds. The mechanism behind this is a first- 
order phase transition induced by pressure. From LDA 
calculations it has been demonstrated [2l[ that in the 
case of zinc-blende compounds with a direct conduction 
band minimum at Ti c and the L\ c state above the Fi c , 
the corresponding wurtzite compound will also have a 
direct gap which is slightly larger. 

In this paper, we present the results of ab initio 
DFT/LDA quasiparticle calculations of band structures 
for both zinc-blende and wurtzite phases of InAs with 
and without the inclusion of the In Ad electrons among 
the valence states. The paper is organized as follows. 



At first we present the technical details of the calcula- 
tions. Then we investigate the role played by the In Ad 
electrons by comparing the US-PP LDA band structures 
calculated with and without the inclusion of the d states 
in the valence band. By means of a further comparison 
with the experimental (zb) energy of the d states we have 
estimated a correction factor for the band structure when 
calculated with the d states frozen in the core. After ob- 
taining the GW band structure in this case, we have ap- 
plied this correction to arrive at the final wurtzite band- 
structure. Finally we treat the case with the d states in 
the valence band by using the screened exchange (SX) 
functional for exchange and correlation Q. This gives a 
band structure with a correct band ordering and hence 
it is a valid starting point for the evaluation of the GW 
corrections. 



III. COMPUTATIONAL DETAILS 

All the calculations are performed with the VASP 
code |22[, a software package for ab initio simulations 
which relies on the description of the electron-ion interac- 
tion via ultra-soft non-norm-conserving pseudopotentials 
(US-PP) f23]or the projector-augmented wave (PAW) 
method [24 . |25[ and a plane- wave expansion of the eigen- 
functions. Both kinds (US and PAW) of pseudopotentials 
are used throughout this study, for both polymorphs and 
for both treatments of the In Ad states. Hence, the num- 
ber of studied cases adds up to 8, for each of which an ac- 
curate convergence study was performed, in such a way as 
to achieve an absolute convergence of a few meV (always 
better than 10 meV) for every electronic level. In this 
paper we use a simplification for the self-energy based on 
a model dielectric function which, in turn, is based on a 
plasmon pole approximation for the dynamics [26l |27| . 

We have performed the convergence study for the two 
polymorphs with the In Ad electrons as either valence (d- 
val) or core (d-core) states. For the zinc-blende phase we 
have used a T-centered Monkhorst-Pack [28[ grid with a 
llxllxll mesh in reciprocal space. For the wurtzite 
phase, instead, a 8 x 8 x 8 mesh was used. 

As kinetic energy cutoffs we used 262 eV when the 
d states are in the valence band (for both polymorphs 
and both pseudopotentials). When the d electrons are 
included in the core we used a cutoff 222 eV and 202 eV 
for the zb and wz cases, respectively. 

The structure optimization has been performed by 
minimization of the total energy In the wurtzite case 
we have assumed the ideal values for the c/a ratio and 
the internal parameter («). The resulting energy versus 
volume curves, normalized to one InAs pair, are displayed 
in Fig. [TJ Here, wc can see that the zb is correctly pre- 
dicted to be the stable phase, since its binding energy is 
more negative than the wz by ~ 18 meV. 

The results obtained from the structure optimization 
of InAs for both zinc-blende and wurtzite phases are col- 
lected in Table [T] together with experimental results for 
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FIG. 1: The normalized total energy versus volume for one 
InAs pair for the zinc-blende (diamonds) and the wurtzite 
(circles) polymorphs. Total energies are calculated by using 
US pseudopotentials. 



lattice constants. The experimental values for wurtzite 
are taken from TEM measurements [l|, while, for the 
zinc-blende, the experimental parameter measured at 
room temperature [29[ has been extrapolated to zero de- 
grees Kelvin. From this comparison we see that the cal- 
culated lattice constants when the d electrons are treated 
as valence states are very close to experimental results. 
Hence, these values have been used throughout our fur- 
ther investigation. 

The results reported in literature for the wz poly- 
morph are obtained using HGH pseudopotentials in the 
d-core approximation [13|: a = 4.192 A, u = 0.3755, 
c = 6.844 A. In the same article and with the same ap- 
proach is reported a = 5.921 A for the zb. More results 
obtained at the LDA level with pseudopotentials (d-core) 
are reported for the zb phase as, for instance, a = 5.04 A 
HI, a = 5.902 A [HI, a = 5.95 A M, a = 5.906 A SI. 



TABLE I: Lattice constants (A) for InAs in the zinc-blende 
and wurtzite phases. US denotes Ultra Soft pseudopotentials 
with the In Ad electrons included among the valence (d-val) 
or core (d-core) states. 





a z b 






experimental 


6.0542 A 


4.2839 A 


6.9954 A 


US d-val 


6.0329 A 


4.2663 A 


6.9669 A 


US d-core 


5.8023 A 


4.1060 A 


6.7051 A 



IV. BAND STRUCTURES 
A. Kohn-Sham and quasiparticle band structures. 

The LDA band structures along high-symmetry lines 
in the Brillouin-Zone (BZ) were obtained by solving the 
Kohn-Sham equation [f| at the equilibrium lattice pa- 
rameters and using US-PP. For the wz-InAs, these are 
shown as continuous lines in Fig. [2] and [3] in the d-val and 
d-core case, respectively. When the In Ad electrons are 
included in the valence band, the calculated band struc- 
ture reveals problems for both polymorphs: a negative 
sp gap is obtained, resulting in a "wrong band ordering" . 
The s-likc r lc state - forming normally the conduction 
band minimum (CBM) - lies below the p-like T^ v (or 
Tg,, and Ti v ) state - constituting the valence band maxi- 
mum (VBM) in other zinc-blende (or wurtzite) semicon- 
ductors. We find Tic — Tu v = —0.346 eV for the zinc- 
blende and Tic — r 6l , = —0.317 eV for the wurtzite energy 
gaps. The only difference between the two InAs phases 
is that the wz exhibits crystal field splitting between the 
r 6 „ and Ti v valence band levels. The calculated value is 
A cr = e(T 6v ) - e(Ti„) = 81 mcV. 

The "negative gap" of the InAs in the zb phase has 
already been observed in other LDA calculations that 
include the d states in the valence band as, for instance, 
reported in Ref. [34[. The main cause of this problem 
is the known limitation of the DFT-LDA to correctly 
describe excited states properties like the band structure. 
Actually, the latter is properly described within many- 
body perturbation theory, which requires a zero:th order 
Hamiltonian for the perturbation series. This is usually 
chosen to be the LDA eigenspectrum. Unfortunately, in 
the present case, the wrong starting Hamiltonian, will 
lead to unreliable GW corrections which actually also 
results in an erroneous band structure as can be seen 
from Fig. d 

When the experimental band-gap is as small as it is in 
the present case (E gap = 0.415 eV at K j35j) sources of 
error other than the fundamental one become important 
like the LDA itself. The LDA becomes less accurate in 
more inhomogeneous materials such as the III-IV semi- 
conductors and even worse in the II- VI semiconductors 
with a higher ionicity. Besides, as it will be discussed 
in the next section, the underestimation of the binding 
enery of the In Ad electrons causes an overestimation of 
the pd-repulsion that contributes to the reduction of the 
band gap. 

For the GW corrections we use a model dielectric func- 
tion for which only the high frequency dielectric con- 
stant Coo has to be specified. For this parameter wc 
have, for both polymorphs, used the experimental val- 
ues eoo = 12.3 (InAs) and = 10.89 (GaAs) [H] cor- 
responding to the zb phase. The resulting quasiparti- 
cle band structures for wz-InAs are presented in Fig. [2] 
and [3] (dashed lines) in the d-val and d-core case, re- 
spectively. We cannot use calculated dielectric constants 
because they are dependent on the LDA gap and, hence, 
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FIG. 2: Quasiparticle (solid lines) and Kohn-Sham (dashed 
lines) band structures of InAs wurtzite with the In 4 d included 
in the valence states: at F a "negative gap" is predicted. 



they are too large for both InAs and GaAs. In the InAs 
case, indeed, the Tic — T±5 V LDA gap is negative and, 
in the GaAs case, it is too small. Anyway, the dielectric 
constants do not differ substantially for the zinc-blende 
and the wurtzite phase, and this allow us to use the zb 
experimental value for both phases. The average value 
for wurtzite is, indeed, close to the zinc-blende value. 
This has been recently demonstrated for InN [l3[ and 
ZnO The almost equal LDA gaps for zinc-blende 

and wurtzite InAs indicate a similar behavior. 



trons. Consequently, the d levels are calculated to be too 
close to the top of the valence band and the pd coupling 
is overestimated within the LDA. This results in a shift 
of the p states towards higher energies and, hence, to 
the gap reduction. The undcrbinding of the d states can 
be estimated by comparing the experimental (zb) and 
calculated positions of the d bands with respect to the 
VBM. The values calculated within the LDA arc 14.3 eV 
and 14.5 eV for the zb and wz structures, respectively, 
compared to an experimental result of 16.8 eV [37j (zb). 
Hence we obtain a scaling factor of 14.3/16.8 = 0.851 
- 85% for the zb and 14.5/16.8 = 0.857 - 85% for the 
wurtzite. We will use this value of scaling factor to reduce 
the overestimation of the pd-repulsion calculated within 
the LDA. 

In order to highlight the role played by the d states 
in the band structure calculations of these materials, we 
have studied both polymorphs with a (US) pseudopo- 
tcntial obtained by assuming the In Ad states to be core 
states while accounting for their influence by including 
self-interaction corrections (SICs) in the underlying cal- 
culations for the free atom (l6| . Using this pseudopo- 
tential and the equilibrium lattice constant obtained in 
the rf-val case, we have calculated the Kohn-Sham band 
structure and found that the correct band ordering is 
re-established. The semiconducting nature of InAs is 
now correctly predicted, with energy gaps being 0.068 eV 
and 0.113 eV for the zinc-blende and the wurtzite phase, 
respectively. Moreover, this calculations forms a good 
starting point for perturbativcly adding GW corrections 
to the band structure in the usual way [H, This 
is what we have done, i.e. we have not attempted self- 
consistencv with resDect to the auashoarticle wave func- 

□ 



B. The pd repulsion and the d-core case. 

The calculation of the band structure for the InAs com- 
pound is complicated by the presence of the d states, as 
discussed previously for systems like Ge, II- VI semicon- 
ductors and some III-V semiconductors [HI HH, H3, 
Hil . [l5| . Indeed, even tough the d states are well sepa- 
rated in energy from the lower lying s and p states with 
the same principal quantum number their wavefunctions 
have considerable spatial overlap which leads to a large 
exchange coupling between these states [{| . A proper de- 
scription requires a pseudopotential which includes the s 
and p states in addition to the d valence states. Such a 
description leads to expensive GW calculations. Another 
possibility is to use a pseudopotential which excludes the 
d states from the valence band. 

Pseudopotential calculations based on the LDA and 
treating the In Ad as valence states give smaller values 
of the gap due to the pd hybridizations of the upper va- 
lence band VBM (p-like states) with the d states. The 
strong pd interaction is mainly due to the fact that the 
LDA underestimates the binding energy of the d elec- 




FIG. 3: Quasiparticle (solid lines) and Kohn-Sham (dashed 
lines) band structures of wurtzite InAs with the In 4 d included 
among the core states. 
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The band structure obtained with the d electrons 
treated as core states misses the effects of the pd repul- 
sion, i.e., the calculated GW band gap is larger than the 
"true" one by the energy corresponding to the "true" pd 
repulsion itself. Hence, we have here estimated this cor- 
rection factor following the method outlined in [l3| in the 
case of InN. By comparing the LDA band structures in 
the d-core and <i-val approaches, it is possible to estimate 
the LDA value for the pd repulsion (A p( j). This can be 
identified with the shift of the p states toward higher en- 
ergies when passing from the former to the latter case. 
Indeed, when the d states are "missing" (i.e. very far 
below the VBM), there is no pd repulsion and it is possi- 
ble to assume that the position of the CBM (s-like state) 
remains unchanged in the two situations. This is because 
in the InAs compound the s-d repulsion is negligible [io| . 
Referring to the schematic representation in Fig. [4] and 
taking the Ai and A2, as defined in this figure, from 
the calculated band structures, we obtain an estimate of 
A p d for both polymorphs, that is 0.41 eV for the zb and 
0.44 eV for the wz phase. 



J-val 



d-core 



p-states (VBM) 
s-state (CBM) - 
p-states (VBM) E 



pd 



FIG. 4: Scheme for the estimation of the LDA value of the 
pd-repulsion (A pd ). 

The pd repulsion obtained in this way, i.e. from the 
LDA band structures, is overestimated with respect to 
the "true" one. To arrive at the "true" value A p( j must 
be reduced by the factor that relates the experimental 
energy difference between the VBM and the d levels to 
the same quantity within the LDA, i.e. by 85% for both 
polymorphs. Hence the true pd-repulsion can be esti- 
mated as A pd x 0.85, leading to 0.35 eV for the zb and 
0.38 eV for the wz. By applying this correction to the 
GW band structure calculated in the d-core approach we 
find that the quasiparticle energy gaps of 0.67 eV (26) 
and 0.79 eV (wz) are reduced to 0.38 eV and 0.41 eV, 
respectively. 



V. LDA-SX-GW APPROACH TO THE BAND 
STRUCTURE CALCULATION 

A. The InAs case. 

In the previous section we have obtained band gaps 
for the two polymorphs of InAs using a method which 
is based on physical intuition but has some theoretical 
limitations. It is based on experimental data for the ri- 
band positions and can thus not be considered as fully 



ab initio. The method does not properly treat the In 
Ad states as valence states as one would require from a 
proper theory. The advantage of the approach, however, 
is the remarkably low numerical effort and the fact that it 
succeeds in reducing the self-interaction among the very 
localized Ad electrons of In. 

We have, therefore, also applied a different approach 
to the problem, which is better defined from a theoret- 
ical point of view. The leading idea is to search for a 
starting point for the GW calculation which is closer to 
the quasiparticle electronic structure, thus approaching 
a self-consistent GW calculation. As mentioned previ- 
ously, in the d-val case the LDA band structure is a bad 
starting point for calculating the GW corrections. More- 
over, even if calculated correctly, it might be inadequate 
to add them in a perturbative way. 

A full quasiparticle calculation in the GW approxima- 
tion involves several difficulties such as the calculation 
of the frequency-dependent dielectric function and the 
solution of a Schrodingcr-likc equation with a non-local 
and energy dependent self-energy operator which is non- 
hcrmitian as opposed to the case of most eigenvalue prob- 
lems involving a local energy independent potential. On 
top of these difficulties the calculations should actually 
be carried out self-consistcntly meaning that the output 
quasiparticle wave functions and energies should be used 
in calculating the self-energy operator from which the 
results are obtained. 

Performing a self-consistent GW calculation is a heavy 
undertaking. It is well known that a non-selfconsistent 
calculation, i.e. approximating GW by Go Wo where the 
screening (Wo) is taken from the random phase approx- 
imation (RPA), gives reasonable results. The problem 
is that there is no systematic way of obtainin g su cces- 
sively better approximations to the self-energy [4l| and, 
hence, for the self-consistency procedure. Indeed, when a 
full self-consistent GW calculation is performed, the self- 
consistent screening (W) is an auxiliary quantity with a 
rather unphysical behavior [42|, . 

In practice, one should always start the many-body 
perturbation expansion from a non-interacting Green 
function with one-electron eigenvalues close to quasi- 
particle energies. Then, self-consistency with respect to 
the eigenvalues can be carried out. Self-consistency with 
respect to the wave functions raises a different problem: 
if the self-energy is energy dependent and non-Hcrmitian, 
the quasi-particle amplitudes are not orthogonal and can- 
not form the basis for a non-interacting Green function 
to start from. Hence, self-consistency with respect to 
wavefunctions should not be attempted. Indeed, when 
this was done, as in Ref. [44], the results were worse 
as compared to those obtained when the self-consistency 
was restricted to the eigenvalues. Within the Screened- 
Exchange approximation (SXA), instead, the self-energy 
operator is Hermitian and energy independent. Hence, 
the wave functions are orthogonal and self-consistency 
can safely be performed. 

Already the founder of the GW approach designed ap- 
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proximations to the GWA |7| in order to simplify its 
application to real solids. One such approximation is 
the Coulomb-Hole plus Screened-Exchange approxima- 
tion (COHSEX) which splits the self-energy of the GWA 
in two parts. The first part , th e Coulomb-Hole term, 
was shown already by Hedin [45l ] to have a rather weak 
dispersion although it is far from negligible. The main 
effect of this term is thus a constant shift of the entire 
band structure. Consequently, neglecting this term en- 
tirely and just keeping the Screened-Exchange (SX) term 
constitutes a rather good starting point for later adding 
the full GW correction in a perturbative way. Also, the 
use of the SX approximation is another way of achieving 
the drastic reduction of the self-interaction among the d 
electrons. 

Viewing the exchange-correlation potential of the LDA 
as an approximation to the self-energy of the GWA and 
believing that the full GWA would predict InAs to be a 
semiconductor, the input LDA wave functions and ener- 
gies are very far from the resulting wave functions and en- 
ergies of the GWA. We here subscribe to the idea that the 
self-energy operator of the SXA being non-local but en- 
ergy independent is much closer to that of the full GWA 
as compared to the LDA potential. The eigenvalues of 
the one-electron quasiparticle equation of the SXA are 
close enough to those of the GWA to allow us to obtain 
the full GWA results from first-order perturbation the- 
ory starting from the SXA. This conjecture is strongly 
supported by recent results by some of us [46j]. We thus 
use the wave functions and eigenvalues from the SXA 
calculation as input for the GWA calculation. 

The calculations were performed by using the model 
dielectric function e(q, p, uj — 0) of Ref. [26[. The dielec- 
tric constant eoo of the material and the average electron 
density p are required to build the static dielectric ma- 
trix. Since the LDA dielectric constant is too large due 
to the zero gap problem, we used the experimental value. 
The static dielectric matrix is then extended to finite fre- 
quencies by the generalized plasmon-pole model [38| . 

The one-electron part of our calculation is based on 
pseudopotentials obtained from the projector-augmented 
wave (PAW) method. This method allows for the recon- 
struction of the all-electron wave functions from their 
pseudo counterparts thus making it possible to calculate 
the important matrix elements of the SX operator be- 
tween the valence states [13, 13 > as we ^ as t° evaluate 
the nonlocal core- valence exchange [l7| ■ 

We would like to stress that this is a very important 
advantage of using PAWs as compared to using ordinary 
norm-conserving or ultrasoft pseudopotentials. Unless 
some additional work toward estimating the correct non- 
local core-valence exchange is done such "normal" calcu- 
lations can easily lead to errors in band gaps of the order 
of several tenths of an eV (47[. The effect is more im- 
portant in the case of valence d states and must thus 
be accounted for in narrow-gap materials with strong 
d contributions to states around the gap - as, e.g., in 
InAs. In the present work, the non-local core-valence 



exchange and correlation are approximated by pure ex- 
change rather than by the dynamically screened poten- 
tial of the GWA. We believe, however, that the effect of 
screening is rather small in the core region. 

We find that the eigenenergies of the SXA calculation 
give the correct band ordering for both polymorphs, i.e., 
positive band-gaps of 0.440 cV (zb) and 0.489 eV (wz). 
Consequently, we have a good starting point - the SX 
eigenspectrum - for the subsequent GW calculation and 
it is now possible to apply first-order perturbation the- 
ory. We obtain a quasiparticle spectrum with an energy 
gap of 0.556 eV for zb and 0.611 cV for wz. The crystal 
field splitting for the wz phase obtained in this approxi- 
mation amounts to 99 meV [49| . The quasiparticle band 
structure of the InAs in the wurtzite phase with the In 
4c? included among the valence states is shown in Fig. [5] 




r K H a r ML 



FIG. 5: GW band structure of InAs wurtzite with the In Ad 
included among the valence states. The quasiparticle correc- 
tions are applied on top of the SX calculation which, in turn, 
is started from the LDA eigenspectrum. 

All the calculations are based on pseudopotentials ob- 
tained in the scalar-relativistic approximation. Hence, 
the calculated energy gaps should be compared with the 
experimental values modified as if there were no spin- 
orbit interaction. In the zb case the spin-orbit coupling 
splits the Ti5 state by shifting the T% v states up by +5A0 
and the T-j v state down by — |Ao. Since the zb gap 
at zero Kelvin is 0.415 eV and the spin-orbit splitting 
Ao = 0.38 eV [29], it follows that the experimental en- 
ergy gap of the zinc-blende in a scalar-relativistic world 
would be 0.542 cV. This is the value which should be 
compared to our calculation. In Table HT1 we see that our 
result compares well with experiment and that this ap- 
proach leads to a better agreement than the three step 
procedure of section HVBI 

Encouraged by the good results for the zinc-blende 
structure, we can estimate what the wurtzite energy gap 
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would be in the real world. Indeed the difference between 
the two polymorphs is mainly in the stacking of the atoms 
(ABC sequence for the zb and A B A B for the wz) , that 
is, they differ in the second nearest neighbor of a given 
atom. Hence, it would not be too far fetched to assume 
that the errors introduced by our approximations will be 
similar in the two cases, thereby allowing us to estimate 
the change in the energy gap due to the change in the 
crystal structure as the difference between the calculated 
wz and zb values, i.e. 

E gap {wz) - E gap (zb) = 55 meV. (1) 

By adding this energy to the experimental gap (at K) 
of InAs in the zinc-blende phase we obtain the following 
estimate of the gap (at K) of the wurtzite polymorph 

E gap (wz) = 470 meV. (2) 

Results from low temperature (5 K) photocurrent mea- 
surements on thick (~ 85 nm in diameter) InAs-InAsP- 
InAs nanowires of Ref Q suggest a value of ~ 550 meV 
for the energy gap of InAs in the wurtzite phase, in good 
agreement with our result. Using this value and by as- 
suming that the spin-orbit splitting for the wz polymorph 
is similar to that of the z&, we find the experimental gap 
of InAs in the wz phase without spin-orbit interaction to 
be 0.677 eV. 



poly \E gap 


LDA 


LDA 


GW 


sx 


GW 


EXP 


type 


d-val 


d-core 


d-core 


ci-val 


ci-val 


w.out S-O 


zb 


-0.346 


0.068 


0.38 


0.440 


0.556 


0.542 


wz 


-0.317 


0.113 


0.41 


0.489 


0.611 


0.677 



TABLE II: InAs energy gaps at F point calculated in the dif- 
ferent approximations outlined in the article: (i) DFT-LDA 
with In Ad included among the valence states (ci-val) (ii) DFT- 
LDA with In 4 c! as core states described by SIC pseudopoten- 
tials (d-core) (iii) GW for the ci-core case after the correction 
for the pd-repulsion (iv) SX calculation with PAW pseudopo- 
tential in the d-val case (v) GW gap resulting from the LDA- 
SX-GW method with PAW pseudopotential in the ci-val case 
(vi) experimental gap without spin-orbit interaction. All val- 
ues are given in eV. 



B. The GaAs case. 

In analogy to the InAs case, we use PAW pseudopo- 
tcntials which treat the Ga 3d electrons as valence states. 
Even though GaAs is correctly predicted to be a semi- 
conductor by the LDA, those band gaps are quite small 
compared to experiment, as can be seen from Table IIIII 
When we compute the quasiparticle corrections starting 
from these LDA eigenfuctions, we find an energy gap 
which is still too small - they are of the order of ~ 1 eV 



in the zb case. A similar result is also reported in the 
case of all-electron GW calculations [13] , indicating that 
one needs to go beyond the non-selfconsistent GWA in 
order to better describe the GaAs band structure. 

We have, thus, applied the procedure described in Sec- 
tion IV Al also to GaAs. The resulting energy gaps are 
summarized in Table [TTTI and the full quasiparticle band- 
structure is reported in Fig.[6]in the wz case. From these 
results we note that, at the LDA level of the calculation, 
the gap of the wz polymorph is 50 meV larger than the 
zb case, consistent with Ref. [2l| and the result for InAs. 
When the SX calculation is performed this behavior is 
reversed and the zb gap is larger than the wz. This situ- 
ation again changes after the final GW calculation, with 
the zb gap being smaller by 219 meV. These calculations 
were performed using the same fc-point mesh as for the 
InAs case [49j . 

It is well known that going from a local potential like, 
e.g., the LDA potential to a non-local and energy inde- 
pendent potential like that of the Hartree-Fock approx- 
imation (HFA) has a drastic effect on band gaps. The 
SX approximation is less extreme than the HFA but we 
still expect a substantial effect on the band gap. This is 
what we actually observe in both semiconductors: in the 
GaAs case the SXA results in a zb gap larger then the wz 
one and, in the InAs case, the SXA calculation changes 
the wrong ordering of the bands into the correct one. 

Besides, we note that there is a delicate balance that 
may be changed in the various approximations of ex- 
change and correlation LDA, SX, GW used. GaAs is 
more sensitive than InAs in this respect. Empirically 
one observes that with increasing size difference between 
cation and anion and increasing ionicity of the bonds the 
wurtzite gap is increased with respect to the zinc-blende 
gap. For GaAs, however, the size difference vanishes and 
the ionicity (compared with ZnO, InN) is the smallest. 

Photoluminescence measurements were recently per- 
formed on single GaAs nanowires using the NSOM tech- 
nique These NWs have mixed crystal structure: 
segments having zb and wz crystal structures can be iden- 
tified in the same wire, as reported in [5l|. Since these 
wires have large wz segments, it was possible to take the 
PL emission from that part of the wire finding that it has 
emission energy lower than that of the zb by ~ 50 meV, 
i. e. the measured energy gap of the GaAs in the wz phase 
at low temperature (~ 10 K) was found to be 1.467 eV. 
Since these NWs are typically 50 nm in diameter, we 
can neglect quantum confinement effects. Our GW re- 
sult is in contrast with these measurements, nevertheless 
we would like to suggest more experiments on GaAs in 
the wurtzite phase, especially experiments performed on 
purely wurtzite GaAs NWs. 

To compare our calculated gaps with experimental 
data we strip the experimental results from their spin- 
orbit contribution, using the same procedure as in the 
InAs case. By using the value 1.519 eV and Ao = 0.33 eV 
for the low temperature energy gap and spin-orbit split- 
ting in GaAs (zb) [2!| we find that the energy gap with- 
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out spin-orbit interaction is 1.629 eV. Assuming that the 
GaAs in the wz phase has a spin-orbit split ting similar to 
the zb one and using the results from Ref. [50] , we found 
an experimental wz gap "without" spin-orbit interaction 
of 1.577 eV. 



poly \E gap 


LDA 


sx 


GW 


EXP 


EXP 


type 










w.out S-O 


zb 


0.330 


1.289 


1.133 


1.519 


1.629 


wz 


0.380 


1.172 


1.351 


1.467 


1.577 



TABLE III: GaAs energy gaps at V point calculated according 
to different approximations: LDA (i), SX (ii) and LDA-SX- 
GW (iii). Experimental gap for the wz [5(|] and the zb with 
(iv) and without (v) spin-orbit interaction. All the calculated 
band gaps are obtained using PAW pseudopotentials and in- 
cluding the Ga 3d among the valence states. All the values 
are in eV. 
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FIG. 6: Band structure of GaAs wurtzite with Ga 3d included 
among valence states. The quasiparticle corrections are ap- 
plied on top of the SX calculation which, in turn, is started 
from the LDA eigenspectrum. 



VI. SUMMARY AND CONCLUSIONS 

We have reported on the first study of the quasipar- 
ticle band structure of InAs and GaAs in the wurtzite 
phase. The calculations were performed within an ap- 
proximately self-consistent GW approach. For the pur- 
pose of comparison we have done this in the InAs case 
also by a three-step procedure based on the LDA in which 
the In 4d electrons were frozen and treated as core elec- 
trons. This resulted in a band structure with the correct 
band ordering. We then calculated the GW corrections 
to this results in a perturbative manner and added a cor- 
rection designed to account for the missing pd-rcpulsion. 
The latter physically rather intuitive approach led to re- 
sults not to far from those of the more fundamental ap- 
proach, which in turn was in good agreement with recent 
experimental results. This demonstrates the importance 
of self-consistency in these materials where the LDA re- 
sults are quite off the mark and actually predict InAs to 
be a metal. 

We have found that the InAs wurtzite energy gap is 
larger than the zinc-blende one by ~ 55 meV, leading 
to the theoretical estimate of the quasiparticle gap of 
InAs in the wurtzite phase (~ 0.47 eV), a value in close 
agreement with very recent measurements on InAs based 
nanowires. 

The quasiparticle energy gap of GaAs in the wz phase 
is also larger than that of the zb polytype. This find- 
ing is, however, in disagreement with preliminary photo- 
luminescence experiment performed on GaAs nanowires 
having a wurtzite crystal structure. 
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